PROCEEDINGS 

OF SCIENCE 



o 

(N 
-i— > 
O 

O 

CO 



D s meson spectroscopy 



Oh! 



> 
00 
CN 

d 



Daniel Mohler , R. M. Woloshyn' 

3 TRIUMFA004 Wesbrook Mall, Vancouver, BC V6T2A3, Canada 
E-mail: jnohler@triumf . ca rwww@triumf . ca 



Preliminary results are presented for the spectrum of D s mesons using the 2+1 flavor Clover- 
Wilson configurations made available by the PACS-CS collaboration. For the heavy quark, the 
Fermilab method is employed and we report on the tuning of the charm-quark hopping parameter. 
As our main focus, we present initial results for the spectrum of P-wave states, where previous 
results have been mostly from quenched calculations. As a cross-check, some calculations of the 
charmonium spectrum are also carried out. 
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1. Introduction 

The spectrum of charmed-strange mesons contains a number of well established states, notably 
the "S-wave" states of quantum numbers J p 0~ and 1 , the D s and the D* and the "P-wave" states 
with quantum numbers 0+ (D* (2317)), 1+ (D,i(2460) and 0^(2536)) and 2+ (D* 2 (2573)). In 
addition there are a number of states which have been observed more recently. These are ^ the 
DJ 1 (2710) the D* 7 (2860) the D* y (3040)) and an unconfirmed state previously observed by SELEX 

the D*y(2632). While the D* { (2710) is commonly believed to have quantum numbers J p = l~, 
there are several possibilities for the other states which are not ruled out by experiment. The 
D* ; (2860) has natural parity and is most often identified with a 3~ state, while some still argue 
the possibility of a + identification [||]. The D* 7 (3040) has unnatural parity and is commonly 
interpreted as either a 1 + or a 2~ state. 

Lattice QCD (LQCD) provides the possibility to elucidate the spectrum without resorting to 
model assumptions. To reach this goal, several systematic sources of uncertainty have to be con- 
trolled. This has recently been achieved for the light-quark ground state meson and baryon spec- 
trum Dealing with heavy charm or bottom quarks on the lattice introduces complications of 
its own, and so far a similar precision for the spectrum of D s states has not been attained. 

One common feature of calculations within various models and early LQCD calculations is 
that the ground states in the + and 1 + channels are often found to be quite a bit heavier than the 
experimental resonances. This has sparked speculations about the nature of these states. Molecular 
or tetraquark interpretations have been suggested. Another contribution to these proceedings [||] 
deals with this possibility. 

In the next section our calculational setup and the tuning of the charm quark mass are de- 
scribed. Some first preliminary results from a limited number of configurations follow in section 
|3| As a cross check, the low-lying charmonium spectrum is also calculated. We conclude with a 
short summary and outlook. 
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Table 1: Run parameters for the PACS-CS lattices [0]. With ciw we denote the heavy quark clover term. # 
configs indicates our current preliminary statistics. 

2. Calculational setup 

Dynamical 2+1 flavor Clover- Wilson configurations generated by the PACS-CS collaboration 
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[Q] are used in this work. They span sea-quark pion masses from 702 MeV down to 156 MeV with 
a lattice spacing of 0.0907(13) fm, as determined in Table [j] shows the parameters for the runs. 
The preliminary results presented here are based on a limited number of configurations as indicated 
in the table. 



2.1 Charm quark treatment 

To determine the mass parameter for the heavy charm quark we use the Fermilab method [J7|] 
as employed by the Fermilab-MILC collaboration [[|] for their efforts involving charm and bottom 
quarks. Within this approach, the charm quark hopping parameter K c is tuned to the value where the 
spin averaged kinetic mass (Md s + 3Mo* ) /4 assumes its physical value. In this simplest formulation 
the heavy quark hopping parameter ce = cb = cj$ is set to its tadpole improved value \, where the 
average link uq is determined from the plaquette. The lattice dispersion relation takes the general 
form [||] 

As there are not enough points to constrain such a fit we use the following simplified fit forms: 



1 neglect the term with coefficient W4 and fit Mi, Mi_ and M4. 

2 fit E 2 (p) and neglect the (p 2 ) 2 term arising from the mismatch of Mi, M2 and M4 

E l (p) &Mf + —p 1 jjpif. 



The left-hand side of Fig. [I] shows an example fit obtained from the first method. The resulting 
parameters are given in the caption. While the statistical errors are somewhat larger with the second 
method, both methods lead to reasonable fits which are compatible within statistical errors. The 
right-hand side of Fig.[j] shows our final tuning results. Based on a linear interpolation, we use 
Kc = 0.12752 for the determination of the mass spectrum. To arrive at this value the results have 
been shifted to account for the somewhat unphysical strange quark mass used in the simulation. 



2.2 Simulation details & Source construction 

Propagators are calculated for 8 source time slices on each (independent) gauge configuration. 
The source locations are chosen randomly within the time slice. For the calculation of the strange 



quark propagators the dfl_sap_gcr inverter from Liischers DDHMC package [|9|, [10p is employed. 
For the charm quark propagators the corresponding inverter without deflation is used. All error bars 
shown are determined with a single-elimination jackknife procedure and are of a purely statistical 
nature. 



A matrix of interpolators is constructed in each channel and the variational method [11, 12] 



used to isolate the low-lying spectrum. To this end, both Jacobi-smeared [13, h4| Gaussian sources 



3 



D s meson spectroscopy 



Daniel Mohler 














^M D +3Mp,)/4 : 







0.127 0.1275 0.128 0.1285 
K charm 



Figure 1: Left-hand side: Example fit for K c = 0.128 obtained from Method 1 (as described in the text). 
The resulting fit parameters are M x = 0.86334(50), M 2 = 0.9337(73) and M 4 = 0.863(28) leading to ^ = 
1.0815(86). Right-hand side: Linear interpolation to determine the K c for the spectrum calculation. Results 
for both the D s 0~ state and the spin averaged mass are shown. The red curve (squares) is a shift to account 
for the effect of a slightly unphysical strange quark mass. 



u s = (S u) x 

N 

S = MS Q with AT = £ K n H n , 

n=0 

H(n,m) = t,( U j (». °) 8 ( H + 1 ™) + U J (" - 7 , 0) + 5 (n - /, m) ) , 

7=1 

and derivative sources Wd, 

D i (x,y)=Ui(x,0)8(x+ly)-U i (x-lO) t 8(x-i,y) , 
W di =DjS 

are calculated. While we use both Gaussian and derivative sources for the charm quark, only 
Gaussian sources are used for the strange quark. 

3. Results for D s mesons and charmonium 

In the 1 + sector interpolators corresponding to different charge conjugation in the mass de- 
generate case are expected to mix. This important effect is illustrated in Fig. ||. While the single 
correlators (left-hand side of the figure) are almost degenerate, one can clearly isolate two distinct 
low-lying states when considering a matrix of two or more suitable interpolators (right-hand side 
of the figure). A similar effect is expected in the 2" channel. 

Figure ^ shows preliminary results for the ground states in the + and 2 + channels. In this 
figure as in all following figures mass splittings relative to the spin-averaged pole mass are plotted. 
The results highlight the need for simulations at light sea quark masses 1 . While our results in 

1 Note that the detailed behavior when approaching the chiral limit is intimately related to the scale setting on the 
lattice. 
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Effective masses from single correlators 



Effective masses from Eigenvalues 
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Figure 2: Left-hand side: Diagonal correlators in the 1 + channel from the K tuning run at K c = 0.128. 
The interpolators corresponding to different charge conjugation in the mass degenerate case lead to almost 
identical masses. Right-hand side: Using a 2 x 2 matrix of correlartors two separate low-lying states can be 
identified. The mixing between the two structures enhances the observed mass splitting. 
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Figure 3: Chiral behavior in the + and 2 + channels. For the latter the results from the T2 lattice irreducible 
representation are shown. The results from the E representation are similar. 



the 2 + channel show only a small deviation from the experimental value (fully consistent with 
the presence of non-negligible discretization effects), the results in the + channel still suggest a 
ground state energy somewhat larger than the physical one. This could be attributed to volume 
effects, discretization effects and/or the lack of scattering states in our basis. 

To get an idea of the magnitude of discretization effects it is instructive to take a closer look 
at the ground state hyperfine splitting. Figure |] shows the results for both D s mesons and for 
charmonium. As expected [||] within the Fermilab approach the discrepancy is smaller for D s 
than for charmonium, where the experimental value is underestimated by about 1 8 MeV. One may 
expect discretization effects for the P-wave states to be of similar magnitude. 

While we have not accumulated enough data points for a sensible extrapolation, one can get a 
first idea of the low-lying charmonium spectrum by comparing our simulation results at the lightest 
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Figure 4: Hyperfine splitting for D s and charmonium compared to the physical splitting. The leading 
discrepancy can be attributed to discretization errors. 
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Figure 5: Overview of low-lying charmonium states. The data points are from the simulation at m K = 156 
MeV and the errors are statistical only. Where appropriate, the error of the experimental state is indicated 
by the thickness of the lines. 



pion mass with the experimental spectrum. Such a comparison is shown in Fig. |5[ 

4. Summary and outlook 

Preliminary results were presented for a simulation of the charmed-strange D s mesons and of 
the low-lying charmonium spectrum. As a tool of choice the variational method has been used. 
The calculations indicate the importance of light sea quarks for some of the states. Results for the 
hyperfine splitting and the low-lying charmonium spectrum show that a reasonable mass-tuning for 
the charm quark has been achieved. For the D s mesons with quantum numbers J p = 1 + the im- 
portance of mixing between different interpolator structures is highlighted. While almost physical 
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sea-quark masses clearly improve agreement with experimental data, D s results in the 1 + and + 
channels, where multiparticle thresholds are close to the state of interest, still show a substantial 
deviation from experiment. The role played by discretization effects, by finite volume effects and 
by an incomplete variational basis (providing only poor overlap with multiparticle states) will be a 
topic of future research. 
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